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ABSTRACT 

We performed a series of high-resolution coUisionless N-body simulations designed to study the sub- 
structure of Milky Way-size galactic halos (host halos) and the density profiles of halos in a warm dark 
matter (WDM) scenario with a non- vanishing cosmological constant. The virial masses of the host halos 
range from 3.5 x 10^^ h~^MQ to 1.7 x 10^^ h~^MQ and they have more than 10^ particles each. A key 
feature of the WDM power spectrum is the free-streaming length Rf^wDM which fixes an additional 
parameter for the model of structure formation. We analyze the substructure of host halos using three 
Rf,WDM values: 0.2, 0.1, and 0.05 Mpc and compare results to the predictions of the cold dark matter 
(CDM) model. We find that guest halos (satellites) do form in the WDM scenario but are more easily 
destroyed by dynamical friction and tidal disruption than their counterparts in a CDM model. The 
small number of guest halos that we find in the WDM models with respect to the CDM one is the result 
of a lower guest halo accretion and a higher satellite destruction rate. These two phenomena operate 
almost with the same intensity in delivering a reduced number of guest halos at z = 0. For the model 
with Rf,wDM =0.1 Mpc the number of accreted small halos is a factor 2.5 below that of the CDM 
model while the fraction of destroyed satellites is almost twice larger than that of the CDM model. The 
larger the Rf^wDM value the greater the size of these two effects and the smaller the abundance of 
satellites. Under the assumption that each guest halo hosts a luminous galaxy, we find that the observed 
circular velocity function of satellites around the Milky Way and Andromeda is well described by the 
Rf,WDM — 0.1 Mpc WDM model. In the Rf^wDM — 0.1 ~ 0.2 Mpc models, the surviving guest halos 
at z = — whose masses are in the range Mh ~ 10^ — 10^^ H'^Mq — have an average concentration 
parameter Ci/5 (= r{Mh)/r{Mh/5)) which is approximately twice smaller than that of the corresponding 
CDM guest halos. This difference, very likely, produces the higher satellite destruction rate found in the 
WDM models. The density profile of host halos is well described by the NFW fit whereas guest halos 
show a wide variety of density profiles. A tendency to form shallow cores is not evident; the profiles, 
however, are limited by a poor mass resolution in the innermost regions were shallow cores could be 
expected. 

Subject headings: cosmology:theory — cosmology: dark matter — galaxies: formation — galaxies: halos 
— methods: numerical 



1. INTRODUCTION 

Non-baryonic dark matter is an essential ingredient of 
current inflation-inspired models of cosmic structure for- 
mation in the universe. From the point of view of par- 
ticle physics, there is no obvious prefe rence for any of 



the predicted dark m atter candidates (Colombi, Dodel 
on, fc Widrow 100(:), which, according to their rmo v e 



seem to be in conflict with observations. 

One of the potential problems of the CDM scenario 
is that the predicted number of low-mass halos — where 
probably dwarf galaxies form — within a Milky Way-size 
halo, greatly exceeds the observed abundance of satellite 
galaxies i n the Local Group (K lypin et al. 1999, hereafter 



locity at the time of their decoupling , can be cold, warm , 
or hot. From the point of view of structure formation, the 
most compelling candidate has been the cold dark matter. 
The CDM acenario for atrueture formation has aucceaafuUy 
Qooountod for oovoral oboorwational footo, partioularly on 



|KKVP99| ; [Moore et al. 1999at see also [Kauffmann et al 
1993). A second problem is that the predicted inner den- 



large scales, without introducing an additional free param - 



sity profiles of CDM halos may disagree with the shallow 
profiles inferred from the r otation curves of dwarf and low 
surfa.ce brightness galaxies (|Moore 1994 ; Florcs fc Primac" 
1994 iBurkert 1995| ; |de"Blok fc McGaugh 1997| pernande: , 
fc Gilmore 199S JTalthough the observational data for the 



eter related to its particle distribution function in phase 
apace. However, on amall acalca and/or in high density re 



latter galaxies are controversial (van den Bosch et al. 1999 
Swaters 



gions of the universe, the predictions of the CDM models 



al. 2000 



Madore, fc Trewhella 200q ; but see Firmani eq 
j^. High-resolution gravitational lensing maps of 
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a cluster of galaxies have also reveale d a soft inner mass 



distri bution in the ha lo of this cluster (Tyson, Kochanski, 



& Del r Antonio 1998 ). The rotation curve decompositions 
of normal galaxies and the TuUy-Fisher relation obtained 
in galaxy formation models as well as the dark mass con- 
tained within the solar radius in our Galaxy, also point 
out to dark halos shallower and/or much less concentrated 
than those predicted by the CDM model lAvila-Reese 
Firme lni, fc Hernandez 1998|; [Navarro 19981 ; |Navarro~ 



Stein metz 1999| ; [Firmani fc Avila-Reese 200C| ; |Mo fc Mac 
20000 If these shortcomings of the CDM scenario are con 
firmed with more observational and theoretical data, new 
alternatives (cosmological and/or astrophysical) have to 
be explored in order to modify the properties of the mass 
distribution at small scales. 

In a recent burst of papers, explored alternatives in- 
clude modifications to: either the nature of the dark mat- 



ter ca ndidate (e.g., Spergel fc Stcinhardt 2000; Hannestad 



1999; gommer-Larsen fc Dolgov 2000|; [White fc Croft 2000 



Firmani et al 2000a ; Hogan fc Dalcanton 2000|; Moore 



et all 2000| ; [Yoshida et al. 200q ; [Burkert 2000 
bles 2p00t Hannestad fc Scherrer 2000; Riotto & Tkachev 
2000), or the generation of the primor dial power spec- 
trum (e.g 
five astrop 



Kamionkowski & Liddle 199£). More conserva- 
lysical mechanisms to overcome the pro blems 
ment i oned above have also been proposed (e.g., [N; 
Eke, i 'z Frenk 1996 ; Gelato fc Sommcr-Larsen 19991; IBuI 



lock, : Cravtsov, fc Weinberg 20'0C[ rp inney, Gerhard, fc Silkj 
2000^ One possible modification is to go from a CDM 
scenario to a warm dark matter (WDM) one. The WDM 
particles (warmons) would suppress the power at small 
scales by free-streaming out of overdense regions limit- 
ing the formation of substructure at scales below the free- 
streaming scale. At large scales, the structure formation 
would proceed in a very similar way to that of a CDM 
model. N-body simulations have shown that indeed large- 
scale structure in WDM models looks similar to that of a 
CDM model (Colombi et al. 1996). On the other hand, as 
Hogan fc Dalcanton (2000) noted, the finite phase density 
of dark halos inferred from observations could be pointing 
to a non-negligible DM velocity dispersion at the time of 
structure formation. 

Using the Press-Schechter formalism, Kamionkowski fc 
Liddle (1999) have shown that if the CDM power spectrum 
is filtered at scales corresponding to dwarf galaxies, then 
the abund ance of Milky Way sa tellites can be reproduced. 
Recently, White fc Croft 2000| reported results from N- 
body simulations for WDM models at high redshifts. They 
found that the abundance of lO^°ft-~^M0 halos is reduced 
by a factor of 5 at z = 3 with respect to the CDM model 
when the power spectrum is filtered at fc « 2h Mpc^^. At 
the same time they showed that the Ly-a power spectrum 
at this redshift is very similar to that of the CDM model, 
which is in agreement with observations. This apparent 
contradictory result is explained by the fact that the col- 
lapse of large-scale structures, as they go non-linear, re- 
generates the init ially suppressed smal l-scale modes in the 
power spectrum (White fc Croft 2000). 

These results encourage us to explore in more detail the 
predictions of WDM N-body simulations at the present 
epoch. Docs the suppression of power at small scales of 
a WDM model actually eliminate the excessive degree of 
substructure predicted by the CDM scenario? Are the 



WDM halos less concentrated? And if so, do they have 
a smoother inner mass distribution than their counterpart 
CDM halos? The main aim of this paper is to give a quan- 
titative answer to the first question. To this end we have 
carried out high-resolution N-body simulations of Milky 
Way-size galactic halos in three different WDM models. 
A host halo of about 2 x 10^^ /I'^M© has more than 10 
particles in the simulations. Since the most successful vari- 
ant of the CDM models is a flat universe with a non-zero 
cosmological constant (JIa = 0.7 and h — 0.7), here we 
also use this cosmological model but instead of CDM we 
introduce WDM with the extra free parameter Rf.wDM- 
these models will be our AWDM models (for economy we 
drop off the greek letter A hereafter when we refer to either 
CDM or WDM models). We will also address the ques- 
tions of concentrations and density profiles of dark halos, 
although the small number of large high-resolved halos and 
the small range of masses in the simulations constrain our 
predictions on this subject. 

In Section 2 we discuss the WDM models to be explored 
in this paper. In Section 3 we briefly describe the numeri- 
cal technique that we used for the simulations. Section 4 is 
devoted to the analysis and comparison with observations 
of the circular velocity function of satellites within host 
halos of Milky Way-sizes. The concentrations and density 
profiles of the host and satellite halos are presented in Sec- 
tion 5. In Section 6 we discuss some of the results, and 
summarize of our main conclusions is given in Section 7. 

2. THE AWDM COSMOLOGICAL MODELS 

Several observational tests such as the distribution of 
galaxies (e.g.. Peacock fc Dodds 199 4), cluster mass es- 
timates (e.g., Carlberg et al. 1996[ ), the determination 
of the baryon fr action in clusters (e.g., [Mohr, Mathiesen 



fc Evrard 1999), and the evolution of cluster abundance 
(e.g.. Eke et al. 1998) point to a cosmological CDM model 
with a low matter density, f2o ~ 0.3. This model also 
successfully accounts for the observationally inferred val- 
ues of the Hubble constant and the age of the universe. 
On the other hand, according to a prediction of the infla- 
tionary theory, the universe should be flat, i.e. a contri- 
bution to the density of the universe from a cosmological 
constant is necessary if JIq ~ 0.3. It has been inferred re- 
cently from observations of high redshift Supernovae that 



the universe is expanding with positive acceleration ( Perl- 
mutter et al. 1999[ ; [Riess 1998[ ^chmidt et al. 1998[ ). Re- 



markably, the estimated value of the cosmological constant 
density parameter is in this case f^A ~ 0.6 — 0.8. Thus, 
the most popular cosmological model has become a flat 
CDM model with a non- vanishing cosmological constant: 
rio ~ 0.3, JIa w 0.7, and h = 0.7 (the Hubble constant in 
units of 100 km s~^ Mpc^^; this value is con sistent with 
the current observational determinations, e.g., Nevalainen 



fc Roos 1998). Here we will use this cosmological model 
but instead of CDM we will introduce WDM. 

A dark matter particle is usually defined as hot or cold 
if at the moment of decoupling from the rest of the cos- 
mic plasma it is relativ istic or non-relativistic, respectively 
( [Kolb fc Turner 1990| ). The classic and only example of 
detected dark matter are the neutrinos. They are hot be- 
cause they were relativistic at the moment of their decou- 
pling. If the mass of the dark matter candidate is much 
higher than 1 GeV and the strength of its interactions is 
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comparable to that of the weak interaction, then it would 
behave as cold dark matter. The thermal velocities of 
these particles at the time of structure formation is neg- 
ligible. In this paper we are interested in a warm DM 
candidate, a thermal relic that at the time of its decou- 
pling was relativistic and whose mass mw is much higher 
than that of its hot counterpart. The Cowsik-McClelland 



bounc prohibits any candidate with a mass larger than 
^ 15 cV (assuming ftwDM ~ 0.3 and h — 0.7) which de- 
couple s when the temperature of t he universe was a few 
MeV dCowsik fc McClelland 1972| ). Thus, the warmon 
should decouple earlier than a hot candidate, in an epoch 
when the total number of degrees of freedom of relativis tic 
particles was certainly very high ( Kolb fc Turner 1990| ). 



Unlike the CDM case, the small-scale density fluctua- 
tions are damped out in a WDM scenario by the free- 
streaming of DM particles. It is straightforward to com- 
pute the comoving free-stream ing scale Rf^wDM (e.g., 
Somnier-Larsen & Dolgov 200C ) : 



R 



•f^wDM =0.2 {QwDMhy/'' (^) '^'mpc. (1) 

The WDM scenario was not attractive in the past be- 
cause of the introduction of the extra free parameter 
Rf,WDM, and because particles in the required mass range 
of ~ 100 eV— 1 keV were not particularly compelling. Nev- 
ertheless, these arguments are somewhat obsolete nowa- 
days. As mentioned in the introduction, the CDM sce- 
nario seems to be in disagreement with observational data 
at small scales, so that models with extra degrees of free- 
dom might be necessary. On the other hand, light WDM 
candidates as palatable as the CDM ones are also predicted 
by particle physics beyond the standard model; one possi 



ble ex ample are the right-h anded neutrinos (e.g., Colombi 
Dodefeon, fc Widrow 1996|) . 



We have represented t he WDM power spec trum by 
the following expression (Bardeen et al. 1986| ; see also 
bommcr- Larsen fc Dolgov 200q ) 



Pi 



WDM 



(fc) 



TwDM{k)PcDM(k), 



(2) 



where the WDM transfer function is approximated by 



TwDKiik) = exp 



kR 



f.WDM 



(kRf, 



WDM ) 



(3) 



and PcDM is the CDM power spectrum which is in turn 



approximated by the following formula (Klypin & Holtz- 
man ^99^: 

P(fc) = 



Ak 



Since our aim is to study the substructure in Milky Way- 
size halos, we simulate WDM models with three differ- 
ent values of Rf^wDM namely 0.2, 0.1, and 0.05 Mpc, for 
which, according to eq. (1), mw is 605, 1017 and 1711 eV, 
respectively. Particle masses of this order were recently 
proposed with the aim to predic t fewer Milky Way satel- 
lites t han in the CDM scen ario ( Kamionkowski fc Liddk 
1999| ; [White fc Croft 20001 ). For a 1 keV particle and for 
flwDM ~ 0. 3, the rms velocity of the p articles is 2km/s 
at z = 40 ( Hogan fc Dalcanton 20do| ). This velocity is 



(1 - 1.5598fcV2 + 47.986fc + 117.77fe3/2 + 321.92fc2)'''"-'''"'' 

(4) 

This power spectrum was obtained by a direct fit to the 
power spectrum estimated using a Boltzmann code and is 
normalized to ug = 1-0, close to the cluster abundance and 
the 4 yr COBE-DMR normalization. Here erg is the rms 
of mass fluctuations estimated with the top-hat window of 
radius 8/i~^Mpc. 

*The number of mass levels was restricted by the amount of memory 



too small to affect the structure of our simulated host ha- 
los. We thus do not consider in our initial conditions the 
thermal component contribution to the velocities of the 
particles. 

3. NUMERICAL SIMULATIONS 

A set of high-resolution simulations of Milky Way-size 
halos (defined as host halos) have been performed us- 
ing a varian t of the Adaptive Refinement Tree (A RT) N- 
body code (Kravtsov, Klypin, fc Khokhlov 1997). The 
ART code achieves high spatial resolution by refining the 
base uniform grid in all high-density regions with an au- 
tomat ed refinement algor ithm. This version of the ART 
code ( Klypin et al. 200C| ) has the ability to handle parti- 
cles with different mass. This is used to increase the mass 
and spatial resolution in few selected halos. Following, 
we describe the way in which the simulations have been 
carried out. 

First, we set the number of mass levels in the mass hier- 
archy to four in all simulations''^. Particles are eight times 
more massive when they pass from one level to the next 
coarser level. For our selected number of mass levels, the 
mass resolution on the finest level corresponds to a box of 
256"^ particles. The realization of the initial spectrum of 
perturbations is done in such a way that this number of 
particles could be generated in the simulation box. The 
size of the simulation box is defined by the requirement of 
high mass resolution and by the total number of particles 
in the finest level. We set a box size of 15 h^^Mpc on 
a side in all simulations which gives a mass per particle 
TUp — 1.66 X 10^ /i^^Mq on the finest level of mass resolu- 
tion. The omission of wavelenghts larger than 15 h~^Mpc 
may affect the characteristic mass and spatial distribution 
of halos, however, we do not expect it to influence the 
halo internal structure (Frcnk et al. 198S). A host halo of 
about 2 X 10^^ h~^MQ resimulated at high resolution will 
have more than 10^ particles within its virial radius rvir, 
defined as the radius at which the average halo density is 
334 times the background density for our selected cosmol- 
ogy, according to the spherical top-hat model. Once the 
mass hierarchy is fixed, we start by running a low-mass 
resolution (LMR) simulation with 32^ particles in a mesh 
with 256"^ cells from which potential host halos are to be 
identified for their future resimulation and analysis. 

Second, the Boun d Density Maxima (BDM) group find- 
ing algorithm (e.g., Colni et al. 1999| ) is used to locate 
potential host halos in the LMR simulation. The BDM al- 
gorithm finds the positions of local maxima in the density 
field smoothed at the scale of interest and applies physi- 
cally motivated criteria to test whether a group of parti- 
cles is a gravitationally bound halo. The friends of friends 

of the computer where the simulations were performed. 
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(a) 




(b) 




(c) 




(d) 




Fig. 1. — The distribution of dark matter particles inside a sphere of 400 h~^kpc radius (solid circle) for the host halo II in four 
different models: (a) CDM model, and (b), (c), and (d) the WDM models with Rf\wDM ~ 0.2, 0.1, and 0.05 Mpc, respectively. 
To enhance the contrast we have color-coded DM particles on a grey scale according to their local density (a pgplot code kindly 
provided by A. Kravtsov) and removed all DM particles whose density was lower than a certain value. The local density at the 
particle positions was computed using SMOOTH, a publicly available code developed by the HPCC group in the UW Department 
of Astronomy. 



group finding algorithm is then used on these halo popula- 
tion to identify the degree of isolation of each halo (halos 
are also visually located). Only halos which are relatively 
isolated^ were taken as potential host halos because we 
are interested in comparing our numerical results with the 
observed substructure in the Milky Way and Andromeda 



galaxies ( |KKVP99| ; [Moore et al. 1999a| ). 

Third, we trace back all particles within a radius r 
1.5 rvir of each selected host halo to get their lagrangian 
positions at z = 40. This radius is large enough to keep 
the contamination due to the presence of particles from the 
second mass level at the level of ^ 1% in mass. We then 



^There are a couple of groups with a mass ^ 10^^ h ^Mq in our box. Milky Way-size halos belonging to one of these two groups were 
rejected for a subsequent analysis. 
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Table 1 

Structural Parameters of Host Halos 



^f,WDM 


Most nalo name tag 


T/ 


J W host 


M 

^ ^ satellite 


cnfw 




62 


T 
1 


(_iVipC ) 




(kni/s) 


[1 X W a Mq ) 












0.2 


I 


258 


3.3 


6 


10.8 


0.38 


0.32 


0.86 




II 


246 


3.1 


7 


9.4 


0.22 


0.19 


0.87 


0.1 


I 


270 


3.4 


13 


12.0 


0.41 


0.29 


0.76 




II 


255 


3.2 


11 


10.5 


0.31 


0.26 


0.85 




III 


241 


2.1 


4 


10.2 


0.34 


0.14 


0.47 




IV 


196 


1.7 


13 


7.6 


0.41 


0.26 


0.69 


0.05 


I 


271 


3.5 


15 


12.6 


0.45 


0.35 


0.83 




II 


258 


3.3 


22 


11.1 


0.37 


0.31 


0.87 




III 


246 


2.1 


13 


13.9 


0.34 


0.18 


0.58 


0.0 


II 


263 


3.3 


35 


11.0 


0.36 


0.34 


0.96 



regenerate the initial distribution using all particles with 
the four different weights (1, 8, 64, 512 xnip). The farther 
away the particle is from the host halo the more massive 
the particle is. The simulation with Rj.wdm = 0.1 Mpc, 
for example, has 1,722,295 particles distributed as follows: 
1,450,496 particles in the first-mass-level, 218,112 in the 
second-mass-level, and 31,040 and 22,647 in the third- and 
fourth-mass levels, respectively. The initial conditions are 
again evolved using ART with the capability of handling 
particles with different mass. The formal force resolution, 
measured by the size of a cell in the finest refinement grid, 
is 0.45 h~^]s.pc and the number of time steps varies from 
325 to 41600. Accurate results are expected at distances 
four times larger than the formal resolution. 

Fourth, the EDM is used once again now to identify 
satellites (guest halos) orbiting around the center of mass 
of host halos. One of the para.meters of RDM is the num- 



ber o f spheres that are randomly placed on the box to 
locate local maxima "seeds" . Wc made sure not to miss a 
significant fraction of guest halos by using the position of 
every fourteenth first-mass-level particle, a number which 
is much higher than the expected number of halos. For ex- 
ample, for the Rf,wDM =0.1 Mpc model we used about 
10^ seeds. Some guest halos have been polluted by more 
than 5% in mass with particles from the second-mass level. 
These guest halos are usually at the periphery of the host 
halos and thus are more susceptible to being contaminated. 
We keep them in our satellite catalogs because we think 
they would still be there even if we increased the size of the 
high-resolution region. To measure the effect on the num- 
ber of satellites due to the smallness of the high-resolution 
volume an additional simulation for one of the host halos 
— the halo II in the model with Rf^wDM = 0.1 Mpc — 
was run doubling the radius of the high-resolution region 
(~ 3 r„). The number of satellites within a sphere of ra- 
dius 200 h~^kpc centered on the host halo is equal to 11 in 
both simulations. This does not mean that there are not 
any differences at all; for example, a satellite, which was 
close to the center 30 h~^'kv>c) in the simulation with 



the sr paller high-resolution region, disappears in the test 
simulation. However, these differences are negligible as far 
as the cumulative circular velocity function for satellites is 



concerned. 

In Table 1 we present the values of some of the physical 
properties of the host halos re-simulated at high-resolution 
for each WDM model (i.e. for each selected R/^wdm 
value). The name tag of the host halo and its maximum 
circular velocity 



(5) 



where M{< r) is the mass of the halo inside radius r, are 
placed in the second and third column, respectively. The 
virial mass and the number of satellites within a sphere of 
radius 200/i~ikpc, centered on the host halo, are displayed 
in the fourth and fifth column, respectively. All satellites 
with more than 10 bound partic les are counted here. The 
conce ntration parameter cnfw (Navarro, Frenk, & White 
1997) is shown in the sixth column. It is defined here as 
the ratio between Vvir and r^, where rg is the NFW scale 
radius. Because our simulations possess the same seed we 
were able to identify a couple of host halos (I and II) in our 
three WDM models and make an inter-comparison study 
of their structural properties. 

We have also measured the ellipticities of the host halos 
using the tensor of inertia. This is defined as 



(6) 



where the sum is over all DM particles within r^ir, Xi 
[i = 1, 2, 3) are the coordinates of the particle with respect 
to the center of mass of the halo, and r is the distance of 
the particle to the halo center. The ellipticities are then 
given by 



ei = l 



Ai 

As' 



62 



1 



A2 

A3' 



(7) 



where A3 > A2 > Ai are the eigenvalues of the tensor 
of inertia. We evaluate the triaxiality parameter using 
the fo llowing formula (e.g.,Franx, lUingworth, & de Zeeuw 



1990) 



T = 



^3 



A? 



(8) 
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A value of T = 1.0 means that the halo is prolate while a 
value of T = 0.0 represents an oblate halo. The elliptic- 
ities ei, 62, and the triaxial parameter for host halos are 
shown in the last three columns. The cUipticies increase 
as Rf.wDM diminishes. In particular, the increment is by 
almost a factor of two for the host halo II when one goes 
from the WDM model with RfWDM = 0.2 Mpc to the 
CDM model. 

Figure 1 provides a visual example of guest halos found 
in our diffferent models, including a CDM model. As the 
host halo we have selected the one denoted by the roman 
number II in Table 1. This halo has a Vmax ~ 250km/s 
but it varies a little from model to model. The letter (a) 
identifies the halo in the CDM model {Rf^wDM = 0) and 
letters from (b) to (d) represent the halo in WDM models 
from Rf^wDM — 0.2 Mpc to R/^wdm = 0.05 Mpc, re- 
spectively. It is notable the absence of substructure and 
a greater roundness morphology of the host halo in our 
Rf,WDM = 0.2 Mpc model as compared with the CDM 
one, which agrees with the analysis of the host halos ellip- 
ticities of the previous paragraph. 

4. THE CUMULATIVE CIRCULAR VELOCITY FUNCTION 
OF SATELLITES 

The first interesting result that arises from our WDM 
simulations is that despite the power spectrum is sup- 
pressed below the free streaming scale, halos of size close 
to or smaller than R/^wdm are formed. This result is not 
obvious at all. Only high-resolution numerical simulations 
could show whether galactic substructures would form and 
survive in such a scenario with a filtered power spectrum 
at high wavenumbers. The number of satellites per host 
halo increases as Rf,wDM decreases (see Table 1). This 
is in agreement with the numerical results at z = 3 of 



White & Croft 2000. Their cumulative number of halos 
above a certain mass increases as the cut-off wavenumber 
ko increases (the power spectrum suffers a sharp drop at 

The present-day cumulative maximum circular velocity 
satellite functions, N{> Knax), for our four models are 
displayed in Figure 2. These functions were estimated as 
follows: for each R/^wdm value we count the number of 
satellites with V^nax greater than a given value within 200 
/i~^kpc from the center of host halos. This number is then 
divided by the number of host halos for each model and 
the volume of a sphere of 200 h~^kpc radius. Although 
we plot this function down to Knax ~ 10 km/s we are 
probably complete only to Knax ^ 20 km/s (KKVP). As 
expected, the number of satellites is much sm aller than 
the one p redicted by a CDM model (KKVP; Moore ct 



al. 1999a| ; Table 1 this paper). The WDM model that 
seems to reproduce better the observed A^(> T^nax) func- 
tion (taken from KKVP) is that with R/^wdm =0.1 Mpc 
or « 1 keV. The discrepancy seen in Figure 2 between 
simulations and observations for Knax ^ 50 km/s can be 
attributed to the intrinsic dispersion of the mass aggrega- 
tion h istory of host halos (Bullock, Kravtsov, & Weinberg 



number of guest halos in a sphere with proper radius 
200 /i^^kpc centered on the host halo at z = 1, which 
is close to the nominal epoch of formation of host halos, 
and compared it with the number we have at z = 0. We 
did this only for the halo 11^ for which a CDM simulation 
had also been performed. At z = 1 we obtain 29, 27, 20, 
and 10 guest halos for R/^wdm = 0, 0.05, 0.1, 0.2 Mpc, 
respectively while at z = the corresponding numbers 
are 28, 13, 9, and 4. Only halos with T^nax greater than 
20 km/s were chosen and this lower limit on Vmax was 
increased by 20% for halos at z = 1 to take into accoun t 
the evolutionary reduction of 14iax (e.g. , [Colin et al. 1999 ) . 




2000) 



V™, (kms/s) 

Fig. 2. — The cumulative maximum circular velocity Vmax 
function for satellites within 200 /i^^kpc from the center of the 
host halo. Solid lines represent the averaged Vmax function for 
each of our model, from top to bottom as Rj.wdm goes from 
0.0 to 0.2 Mpc. The averaged Vmax function from satellites of 
Milky Way and Andromeda is represented by stars (taken from 
KKVP). Error bars are just Poisson errors. For clarity, we only 
plot error bars for the CDM model and the WDM model with 
Rf.wDM = 0.1 Mpc and the observed Vmax function is dis- 
played with dotted error bars. The square at Vmax ~ 80 km/s 
has no error bar because there is only one guest halo at this 
bin. 

The number of guest halos for the CDM model remains 
approxi mately constant with in this proper volume after 
z = 1 (Moore et al. 1999a). Does this mean that the 
guest halos inside this volume are not destroyed and that 
small halos outside this volume do not fall later into this 
volume? We have done the following ex periment to mea - 
sure the degree of satellite destruction (Kravtsov 2000): 
we have tagged all DM particles within 200 /i^^kpc from 
the center of the host halo at z = and traced back them 
at z = 1. We then apply the EDM algorithm on this sub- 
set of DM particles at z = 1. The difference between the 



Is the reduced number of satellites at present time in 
the WDM scenario caused only by the suppresion of power 
at small scales? To answer this question we counted the 

®The BDM algorithm finds two progenitors of the host halo I of comparable mass at 2: = 1, so we decided not to use this halo for the 
argument developed in the paragraph. 
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number of halos found at z = 1 and its number at z = 
gives us the degree of satellite destruction. The percentage 
of destruction is about 35% for the CDM model while it is 
63% for the WDM model with Rf^wDM = 0.1 Mpc. Since 
there is indeed destruction of guest halos from z = 1 to 
z = 0, the only manner to maintain a constant number of 
satellites in the CDM model is by incorporating small ha- 
los from outside the chosen volume (actually, some small 
halos might also form within the volume; however, the 
probability of such an event is very low). The accretion 
rate of guest halos is less in the WDM models just be- 
cause the number of available halos that could fall into 



average cnfw ~ 12. 

The guest halos (Mf, « 10^ - 10" /i-^M©) present a 
wide diversity of density profiles. Unfortunately, since 
the number of particles in these halos is not very large 
(~ 100 — 200 particles for most of them), the resolution 
is not sufficient to study the inner density profile with 
accuracy (see a discussion on this subject in §6.2). For 
those guest halos (~ 15%) whose density profiles are rea- 
sonably well described by the NEW fit we obtain a mean 
Cnfw ~ 8 and 10 for Rf^wDM = 0.2 and 0.05 Mpc models, 
respectively. We estimate a mean Cnfw of 25 by ex trapo 



lating the results o f the CDM model to low masses ( Avila- 



the volume is lower than in the CDM case. We have also Reese et al. 1999) or by using our own CDM simu 



detected that guest halos are more easily destroyed in the 
WDM models. These two effects which are of comparable 
size work together to deliver a much smaller number of 
satellites at z — 0. 

The more efficient disruption of satellites in the WDM 
scenario is very likely due to the fact that guest halos are 
less concentrated in this scenario (see next section). In 
fact, we have found that the guest halos that survived 
until z = are more concentrated than those at z — 1. 
On average, the guest halos at z w 1 — 1.5 are twice less 
concentrated than those at z = 0. 

5. CONCENTRATIONS AND DENSITY PROFILES 

We define the concentration parameter c^/s as the ratio 
between the halo radius r/j^ and the radius within which 
1/5 of the total halo mass Mh is contained. In Eigure 3 
we plot this parameter versus the halo mass for host (large 
symbols) and guest (small symbols) halos for our WDM 
model with R/^wdm = 0.2 Mpc (triangles) and the CDM 
model (crosses). Only those halos which have more than 
90 particles inside their radii have been analyzed. The 
solid and dashed lines in Eigure 3 are extrapolations to 
small masses of the relations Ci/5 — Alh found by Avila- 
Reese et al. (1999) for the corresponding CDM model for 
isolated halos and halos in groups, respectively. From Fig- 
ure 3 one can see that the concentration of host WDM 
halos is only slightly smaller than that of CDM halos. For 
the small guest halos, the difference is more notable; in 
the 10^ — 10^^ H^^Mq mass range the concentrations for 
the case Rf^wDM = 0.2 Mpc are approximately 1.8 — 1.2 
times smaller than those obtained in the CDM model. It is 
interesting to see that the extrapolations showed in Figure 
3 do actually agree with the concentrations of the guest 
halos in our fiducial CDM model. 

The density profiles of host halos obtained in our sim- 
ulations with Rf^wDM = 0.2, 0.1, 0.05, and 0.0 Mpc 
{Mh « 1 — 3 X 10^^ h^^Mo) are shown in Figure 4 (sym- 
bols). In order to avoid overlapping, the profiles from the 
Rf,WDM = 0.1, 0.05, and 0.0 Mpc models were shifted 
in log/9 by —1, —2, and —3, respectively. These density 
profiles are well described by the NEW formula (lines), al- 
though the inner slope in some cases of the WDM models 
is slightly shallower than r~^. The inner profile for the 
CDM halo is actually steeper than r~^. The correspond- 
ing NFW concentration parameters for the host halos are 
given i n Table 1. For compar ision, in the CDM simula- 
tion of Avila-Reese et al. 199S a 10^^ h^^Mr;) halo has in 



ation. 




10" 
Mass (h^^Mg) 



Fig. 3. — The concentration parameter Ci/5 versus the mass 
for host halos (large symbols) and guest halos with more than 
90 particles (small symbols) in the WDM {Rf^wDM = 0.2, 
solid triangles) and CDM (crosses) simulations. The solid 
and dashed lines correspond to the linear fittings found in a 
CDM simulation for isolated hal os and halos within grou p- and 
galaxy-size halos, respectively (Avila-Reese et al. 1999). The 
only host CDM halo here resulted less concentrated than the 
average. 



6. DISCUSSION 

6.1. The accretion and destruction rates of guest halos 

Our results show that in a WDM model with 
Rf,WDM = 0.1 Mpc the number of guest halos within 
Milky Way-size halos agrees with the observed number of 
satellites in the Milky Way and Andromeda galaxies. Us- 
ing the extended Press-Schechter formalism and assuming 
that the amount of substructure is preserved since the host 
halo formation epoch, Kamionkowski & Liddle (1999) also 
found good agreement with the observations for models 
with a sharp drop in the power spectrum on small scales 
and with zero velocity dispersion (they were interested in 
studying the broken-scale-invariant model). However, the 

'^The halo radius r^j is dcflined as the minimum between rvir and the truncation radius (where the spherically averaged outer density profile 
flattens or even increases). In fact, all the host halos and most of the guest halos (more than 70%) attain their virial radius. 
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cut-off length-scale in the power spectrum they requiered 
is two to three times smaller than the one used here in 
order to fit the observational data. We argue that at the 
basis of this difference is the fact that for models with a 
cut-off in the power spectrum and with negligible particle 
dispersion velocities, the amount of substructure in a host 
halo is not preserved. 




Fig. 4 The density profiles of the host halos in the WDM 
models with Rf.wDM ~ 0.2, 0.1 and 0.05 Mpc, and in the 
CDM model {Rf^wDM = 0.0). The profiles of the three lat- 
ter simulations were shifted in the log of the density by —1, 
—2, and —3 in order to avoid overlapping. The different lines 
are fits to the data using the NFW profile and have the same 
correspondence with the symbols in each one of the models. 
The corresponding NFW concentration parameters are given 
in Table 1. 

In Section 4, we have shown that for the CDM model the 
number of satellites within a sphere of proper radius 200 
h^^kpc centered on the host halo remains approximately 
constant from the nominal ep och of host halo form ation up 
to the present time (see also Moore et al. 1999a). Never- 
theless, as discussed in § 4, this is a particular case in which 
the destruction and accretion/formation rates of guest ha- 
los balance up. In WDM models, the accretion/formation 
rate of guest halos is smaller and the destruction rate is 
more efficient than the one in the CDM scenario. These 
two facts operate at the same time to produce a number 
of guest halos at z = much smaller than the one found 
at the nominal epoch of host halo formation. Therefore, 
Kamionkowski & Liddle (1999) overestimated the filtering 
scale of the power spectrum because they did not consider 
that a high fraction of satellites are destroyed by dynam- 
ical effects during their lifetime in the dense environment 
of host halos. Recently, White & Croft (2000) have re- 
ported from numerical simulations of WDM models that 
at z = 3 the number of small halos is reduced roughly by a 
factor of five when the CDM spectrum is filtered at scales 
much larger than 0.1 Mpc. If one takes into account that 



a significant fraction of these halos will be destroyed by 
z = and not many small halos will be accreted, then the 
agreement with observations will be reached for a cut-off 
scale smaller than the one suggested by these authors. 

Why a significant fraction of guest halos become de- 
stroyed in the WDM scenario? We offer the following 
answer: in a WDM model these halos form later than 
in a CDM model, when the mean density of the uni- 
verse was lower. For example, the formation redshift (de- 
fined as the redshift where (7(Mvir,z) = 1) of a halo of 
Mvir = 2 X IO^/i-^Mq is 4.1 for the CDM model, whereas 
for the WDM model with Rf,wDM = 0.2 Mpc it is 2.5. 
Thus, the characteristic overdensity Sc for guest halos of 
this mass in this WDM model is about 4.5 lower than the 
corresponding for guest halos in the CDM model. If Sc 
approximately scales as the cube of the concentration pa- 
rameter, this would mean that these halos were about 1.7 
less concentrated at their formation epoch. This interpre- 
tation agree with the results from our study on concentra- 
tion parameters (§5). In summary, guest halos in WDM 
models are more easily disrupted because of their puffier 
density distribution 

6.2. Density profiles of the WDM halos 

We have found that the present-day density profiles of 
Milky Way- size galactic halos in a WDM model with a 
Rf.WDM value corresponding to masses of the order of a 
few 10^ M0 do not differ from those obtained in CDM 
models. A similar result was also reported by Moore et 
al. (1999b). They cut the power spectrum on a mass scale 
smaller by more than a factor of ten than the mass of the 
studied cluster. The question about how the ffitering of the 
power spectrum affects the inner density profiles of halos 
with masses close or smaller than the cut-off mass scale is 
still open (see Avila- Reese et al. 1998). Unfortunately, we 
could not explore this question in detail in this paper be- 
cause our guest halos with masses close to 2 x 10^ Mq are 
poorly resolved (a study on this direction is in progress; 
small halos in a WDM cosmology are being resimulated 
with high-resolution). Nevertheless, we have found that 
the small WDM halos are actually less concentrated than 
the CDM ones. The difference in the mean in our extreme 
Rf,WDM = 0.2 Mpc model is less than a factor of two in 
the Ci/5 parameter. We estimate that this factor is too 
small to offer a solution to the problem related to the in- 
ner structure of dwarf galaxy dark halos. Thus, although 
the WDM scenario helps, apparently it does not solve this 
problem. Nevertheless, some questions that remain open 
might change this conclusion. In the following, we discuss 
them. 

In the monolithic collapse the orbital tangential veloc- 
ity of the collapsing particles plays a significant role in the 
final virialized configuration; it is expected t hat in a hot 
monolithic collapse shallow cores form (e.g., |van Albada 



1982| ; [Aguilar fc Merrit 1990| ; [Firmani et al. 2000b| ). To 
produce a hot monolithic collapse in a cosmological sce- 
nario, dark matter particles should have a non-negligible 
velocity dispersion (thermal energy) at the redshift of the 
halo formation — there should also be a cut-off in the 
power spectrum at the scale of interest. On the other hand, 
if this is the case, then the formation of shallow cores with 
a maximum limiting d ensity determined by the velocity 
dispersion is expected ( Hogan fc Dalcanton 2000| ). These 
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authors estimate that shallow cores, which are in agree- 
ment with observational determinations in dwarf galaxies, 
can be produced if warmons have a mass mw ~ 200 eV; 
in this case the particle dispersion velocity is significant. 
Of course, according to our results, a WDM model with 
mw ~ 200 eV would produce too few guest halos to agree 
with observations. However, if the inclusion of gas within 
the guest halos plays an important role to avoid their dis- 
ruption, then a larger fraction of guest halos could survive 
until the present epoch. It is probable that with an mw 
only slightly smaller than IkeV, shallow cores that explain 
the rotation curves of dwarf galaxies form if the velocity 
dispersion is taken into account. Nevertheless, as we have 
shown, in this case shallow cores in the massive (host) ha- 
los will not be produced. Thus, if the increment of the core 
radius with the halo mass or T^nax inferred by Avila-Reese 
et al. (1998) and Firmani et al. (2000b) from observa- 
tions is confirmed, then the coUisionless WDM scenario is 
ruled out. It is interesting to note that in a self-interacting 
WDM model, the particle velocit y dispersion would be 
larger than in the coUisionless case (Hannestad & Scherrer 
2000|] A detailed numerical study of this model, which 
appears to be more palatable than a self-interacting CDM 
model, is desirable. 

Finally, we should comment that the problem of an 
apparently excessive number of guest halos in the CDM 
model holds if one assumes that each guest halo hosts a 
luminous galaxy. This assumption was recently challenged 
by Bullock et al. (2000) who proposed that reionization 
can efficiently inhibit dwarf galaxy formation. 

7. CONCLUSIONS 

Using high-resolution N-body simulations, we have stud- 
ied the substructure inside Milky Way-size halos in a 
WDM cosmological scenario. We have also addressed the 
question whether the density profiles of host and guest ha- 
los are different from their corresponding CDM ones. Our 
main conclusions are: 

1. Despite the fact that the power spectrum of fiuctua- 
tions is suppressed at small scales, a non-negligible number 
of virialized structures, corresponding to these or smaller 
scales, form and survive within larger structures. The ac- 
cretion rate of small halos is found to be less in the WDM 
scenario than in the CDM one. This is simply explained by 
the fact that a smaller number of small halos are avalaible 
for their incorporation into host halos in the WDM models. 
A higher satellite destruction rate is found in the WDM 
scenario as compared with the one in the CDM model: it 
can be accounted for the fact that guest halos are less con- 
centrated by about a factor two in average in the WDM 
models. The less efficient halo accretion and the higher 
satellite destruction have almost the same weight as far as 
the final count of satellites within host halos at z = is 



concerned. The larger the Rj^wdm value the greater the 
size of these two effects, and the smaller the abundance of 
satellites. 

2. The predicted maximum circular velocity function 
of guest halos that seems to best fit the observed one for 
satellites in the Milky Way and Andromeda is that given 
by the Rf,wDM = 0.1 Mpc model. This Rf^wDM value 
corresponds to a warmon of mass about 1 keV. 

3. For the R/^wdm — 0.1 and 0.2 Mpc models, guest 
halos (Mh ~ 10^ — 10^^ h~^MQ) have a concentration 
parameter ci/5 which is roughly twice smaller than that 
of the CDM halos. For those guest halos whose density 
profiles are reasonably well described by a NEW paramet- 
ric fit (~ 15%), the cnfw parameter is roughly 1.5 — 3.0 
times lower than that of the CDM halos. This difference 
in the concentration parameters, for both ci/5 and cnfw, 
vanishes as we go to more massive halos. 

4. The density profile of the host halos (Afvir ~ 1 — 
3 X 10^2 h-^MQ) is weU described by the NEW profile 
(in some cases, the inner slope is slightly shallower than 
r~^). The guest halos have a wide variety of density pro- 
files and those whose masses are below the corresponding 
cut-off mass scale probably present a small shallow core. 
The poor mass resolution of the simulations at these scales 
limits our predictions. 

In summary, we have shown that in the WDM model 
with Rf^wDM ~ O.lMpc or mw ~ 1 keV the degree of 
substructure within a Milky Way-size halo is much lower 
than in the CDM model and is in agreement with observa- 
tions, if one assumes that each guest halo hosts a luminous 
galaxy. The problem of cuspy halos probably still persists 
in the WDM scenario, although we have found that halos 
— in particular the small ones — are less concentrated than 
the corresponding CDM halos. If the inclusion of baryonic 
matter helps to significantly avoid disruption of substruc- 
ture, then the agreement with observation may continue 
for less massive WDM candidates for which the rms veloc- 
ity is larger. In this case, the formation of shallow cores 
in dwarf galaxy like halos is expected. Nevertheless, as we 
have shown, halos of larger masses will not have a shallow 
core. 
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